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Abstract 



Defects which are predominant in a reaUstic model, usually spoil its integrability or solvability. 
We on the other hand show the exact integrability of a known sine-Gordon field model with a defect 
(DSG), at the classical as well as at the quantum level based on the Yang-Baxter equation. We 
find the associated classical and quantum R-matrices and the underlying q-algebraic structures, 
analyzing the exact lattice regularized model. We derive algorithmically all higher conserved quan- 
tities Cn, n = 1, 2, ... of this integrable DSG model, focusing explicitly on the contribution of the 
defect point to each C„. The bridging condition across the defect, defined through the Backlund 
transformation is found to induce creation or annihilation of a soliton by the defect point or its 
preservation with a phase shift. 

PACS:02.30.1k , 11.15.Tk , 02.20.Uw , ll.lO.Lm , 72.10 Fk 

Keywords: Sine-Gordon model with defect; classical and quantum integrability; Yang-Baxter equation; 
infinite conserved quantities; soliton creation/annihilation by the defect. 

1 Introduction 

Systems vi^ith defects and impurities are prevalent in nature. Many theoretical studies are dedicated to 
various models with defects starting from classical and semiclassical to quantum as well as statistical 
models [1], with several of them devoted exclusively to the sine-Gordon (SG) model with defect (DSG) 
[2] or inhomogeneity [3], which have enhanced physical importance [4, 5]. A specific form of DSG model 
with a defect at a single point x = : 
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exhibiting intriguing properties close to integrable systems was investigated in a series of papers 
[2, 6]. These investigations aimed to find out mainly the additional contribution of the defect point 
to the conserved Hamiltonian and momentum of the system, using the Lagrangian formalism and 
the effect of the defect point on the soliton solution using the scattering theory. Though the studies 
were concentrated basically on the classical aspects of this model, some semiclassical and quantum 
arguments were also put forward [6]. The important central idea of this approach is the existence of 
an auto Backlund transformation (BT) frozen at the defect point x = 0, relating two solutions u"^ of 
the SG equation along the positive and negative semi-axis [2, 6]. 

Our aim here is to establish the suspected intcgrability of the above DSG [2], by showing the 
existence of infinite set of its conserved quantities and finding them explicitly. The idea is to adopt 
the monodromy matrix approach expressed through matrix Riccati equation [7], a true signature of 
the integrable systems [8] , and couple it with the important concept of extending the domain of defect 
fields u"^ through BT [9]. This approach yields apart from finding out systematically the defect- 
contribution for all higher conserved quantities, an intriguing possibility of creation or annihilation of 
soliton by the defect point. 

More significantly, exploiting the ancestor model approach of [10] we find an exact lattice regu- 
larized version of the DSG model by using the realizations of the underlying algebra. This allows to 
solve the long awaiting problem of establishing the complete integrability of this model, by finding 
the classical and the quantum R-matrix solutions and showing the exact solvability of the classical 
and quantum Yang-Baxter equations (YBE). The exact algebraic Bethe ansatz solution can also be 
formulated for the quantum DSG model, though its explicit resolution needs further study. 



2 Bridging condition and Lax pair for the SG model with defect 

We focus on the central idea in DSG model (1.1) of gluing its fields across the defect point through 
the BT as 



where 



ut{x,0) =Ut {x,Q)+p{x) + q{x), uf {x,0) = {x,Q) + p{x) - q{x) (2.2) 

p{x) = asm — \^ q = a-^ sin — \^ (2-3) 

with parameter a signifying the intensity of the defect. However, we would like to stress on an 
important conceptual diff'erence in the role of BT (2.2, 2.3) played in the present approach and that in 
the previous studies [2, 6], where the above BT was considered to be frozen at the defect point x = 0, 
and hence playing no role at any other point: x ^ 0. Therefore, since the solutions u^{x) can not be 
related through BT at other points along the axis, soliton number remains unchanged while moving 
across the defect point [6] . We on the other hand implement here the idea of [9] used in the semi-axis 
SG model, where the domain of the field is extended with the application of BT. Therefore in place 
of a frozen BT we use (2.2) effectively at all points of the axis including the defect in the following 
sense. 

Define a solution of the SG equation with rapidly decreasing initial data 

/ n\ i u~{x,0) if x < , . ^. { u^{x,0) if x<0 
u(x,0) = < , ) : , and ut(x,0) = < ' , 
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satisfying the gluing conditions (2.2) and having the hmits hm|j.|^^ u{x, 0) = 0. This field solution 
of the SG equation allows to extend the pair of functions u~{x,0), u^{x,0) smoothly from the left 
half-line x < onto the whole line using the BT (2.2) with a limiting value at the positive infinity 
X +00 : u~{x,0) 27rm_ with an integer m_. Following [9] one can prove also the existence 
and uniqueness of such an extension. Similarly one can prolong the other pair of functions ?i+(.'r,0), 
uf{x, 0) from the right half-line to the whole line by means of the same BT and get u'^{x, 0) 27rn+, 
at a; — — oo, n+ being another integer. Now one has two potentials n'*"(x,0), u'f{x,0) and u~{x,0), 
u^{x,0) related to each other by the BT. If the function u{x,t) satisfies the DSG equation then the 
functions u'^{x, t) and u~{x, t) solve the usual SG equation. However in the context of the DSG which 
is the focus model here, such extensions can be considered to be virtual and used for mathematical 
manipulations, while the physically observable fields are only n~ in the domain x <0 and similarly u'^ 
in X > 0. Therefore any solution u~ moving from the left along the axis x < would be transformed 
after crossing the defect at x = to a solution n"^ in the region x > 0, determined through the relations 
(2.2). Therefore, as we see below, it opens up the possibility of creation or annihilation of soliton by 
the defect point, which was prohibited in earlier studies due to consideration of a frozen BT relation 
[6]. Apart from these solutions, a single soliton suffering a phase shift, while propagating across the 
defect point, as found earlier [2, 6], seems also to be present. Interestingly, the BT expressed through 
scalar relations (2.2) can be incorporated more efficiently into the machinery of integrable systems by 
representing it as a gauge transformation relating the Lax pairs of the DSG: 

U{u+) = F^U{u-){F^)-^ + F^{F^)-\ V{u+) = F^V{u-){F^)-^ + Fj^iF^)-^ (2.4) 

where F°(^, u'^,u ) is the Backlund matrix (BM) 

F0(e,u+,u-) = e-i'^3«"M(C,a)ei-3-+, M{C,a) = "^j, (2.5) 

involving both fields and bridging between them at all points, including the defect point x = 0. 
We can check directly from the matrix BT relations (2.4) that by inserting the explicit form of SG 
Lax operators [8]: 

If u , . u \ 1 f u , ■ u \ 

U = — I «to-3 -I- fci cos -(72 + ko sm -ai I , ^ = ^ \ UxCT3 + ko cos -(72 + sm -ai I , (2.6) 

where A:o = C+ |) = ^ — ^, with spectral parameter S^, and comparing the matrix elements, 
one can derive the scalar BT relations (2.2). It is also obvious from (2.4) using the flatness condition 
Ut — Vx + [U,V] =0 that if u~ is a solution of the SG equation, so is u'^. Note also that since 
the corresponding Jost solutions are related by $(^, u"^) = F^{^,u'^ ,u~)^{^,u~), the exact N-soliton 
solution may change its number by one, after crossing the defect point, a possibility lost for the frozen 
BT [6]. 



3 Conserved quantities for DSG model 

For deriving the infinite set of conserved quantities, an essential property of an integrable system, for 
the SG model with a defect we combine the matrix Riccati equation technique for the standard SG 
model [7] with the idea of bridging scattering matrices through BT [9] . Therefore let us first describe 
briefly the technique developed by Faddeev-Takhtajan for the SG model. 
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3.1 Conserved quantities for SG model 

Define the monodromy or the transition matrix as a solution to the associated hnear equation 

dT 

— {x,y,0 = U{x,OT{x,y,0, (3.7) 

with the initial data T{y,y,(^) = 1. To expand the transition matrix in asymptotic power series as 
1^1 GO, it is convenient to gauge transform the variable 

T ^ f{x, y, e) = ^-\x)T{x, y, O^iy) (3.8) 

in equation (3.7), with ii{x) = e4"(^)'^3 g^j^j^j represent it as a product of amplitude and a phase 

f{x, y, 6 = (1 + Wix, 0) exp(Z(x, y, 0)(1 + W{y, OT', (3.9) 

where W{x, ^) is an off-diagonal and Z{x, y, ^) a diagonal matrix, satisfying the condition Z{x, y, $,)\x=y 
0. By a direct substitution of the gauge transformed T in the form (3.9) into the equation (3.7) one 
gets 

Z{x,y,0 = \. J\di^>^ + (C^2 - ^a2e-(->3)W^(x',0)cia;', (3.10) 
where 9{x) = ut{x,t) + Ux{x,t) and W{x,$,) solves a matrix Riccati equation 

^ = Y^easW + ^^((72 - Wa2W) - ^(^^26^"^^^ - Wa2e''"'^W), (3.11) 

This nonlinear equation due to very special form of the coefficients admits asymptotic integration at 
C ^ oo, 

W{x,0 = ^^=o^^, (3.12) 

where Wq = iai. Putting expansion (3.12) in (3.11) and comparing the coefficients with different 
powers of ^ we get the recurrence relation 

dW (x) i 
Wn+i{x) = 2if^3^^ - e{x)Wn{x) + -T.l^^Wk{x)aiWn+i-k{x) - (3.13) 

~^lZlWk{x)c7^e'<->^Wr,-^^k{x) - ^cTie--35„,i, forn = 0, 1, ... 

The corresponding expansion for Z{x,y,^) = crs + iS^;^ ^"^n'^^ ; yields from (3.10): 

Znix, y) = \f^ a2(e^«(^')'^3^„_i(x') - Wn+i{x')dx', (3.14) 

where matrices Wn, Z^ are of the form 

Wn{x) = -Wn{x)a+ + Wn{x)a-, Zn{x) = ^{Znix){I + (73) + -Zn(x)(7 - (73)) (3.15) 

and relations above can be written as the recursion relations starting from n> 1 

Wn+lix) = - 0{x)Wn{x) + ^T.l^^Wkix)Wn+l-k{x) - (3.16) 

-|5]^=o^fc(^)e'"(")^n-i-fe(^) - |e-(5„,i, 
and Zn{x,y) = - f {wn+i{x') - e-'^'^'^ Wn-i{x'))dx' . (3.17) 

4 Jy 



with wq = i. To derive finally the set of conserved quantities we take the limit of the monodromy 
matrix f{x,y,Ox-^+^,y^^oo = TiO = e^(«) + 0(|^|-°°) where 

m = ^ipm + ^3) - pm - ^3)), p{o = ^_^^iun_^_j^^=,^^^ - ^(x - y)). (3.18) 

As shown in [7] the generating function of the conserved quantities : p{^) = loga(^) = iSJ^j^^, at 
1^1 — > oo. is obtained by solving the recurrence equation (3.16) as 

Ci = --/ {-{utix)+ux{x)f + {1- cos u{x))dx (3.19) 

4 J —oo ^ 

and for arbitrary n > 1 

Cn = - {wn+i{x) - e-*'^(^)«;„_i(a;))da;. (3.20) 

4 J— oo 

To derive the asymptotic expansion for ^ ^ it suffices to use the involution [7] tt, u) (— tt, —u), 
with TT = Ut, which leaves the Lax pair invariant. As a result we get loga(^) = iSJ^j^C-n^", as ^ ^ 
0, where Cq = ^ lim^i^+oo ^^(^^5 i) and C_,„(-7r,?i) = (— l)"C„(7r, — u), n = 1,2,..., giving in particular 
C_i = — J Jl^^ {^(utix) — Uxix))"^ + (1 — cosu{x))dx. Therefore one can get the explicit form of the 
momentum P and the Hamiltonian H of the SG model as 

/oo 
P{u)dx, P{U) = UxUt, 
-oo 

/oo 1 
H{u)dx, H{u) = -^{ul + u^) + (1 - cos-u). (3.21) 
-oo ^ 

3.2 Extension to DSG model 

We now extend the above result of the standard SG model to the SG with a defect (DSG) showing 
that the DSG equation admits an infinite set of conserved quantities indicating the integrability of this 
system. In fact, any conserved quantity C„ = Pn{x,t)dx of the SG model can be transformed 
into a conserved quantity for the DSG model by adding some extra term D„, as the contribution from 
the defect, such that 

/O r+oo 
Pnix,t)dx + D„+ / pn{x,t)dx. (3.22) 
-oo Jo 

Our aim is to find an algorithm for evaluating the additional terms Dn, for which we suitably modify 

the above approach for the SG model [7] by using (2.4), a crucial relation in the DSG model. 

In analogy with the SG we define the monodromy matrix of the DSG as a solution to the associated 
linear equation with a defect at the point x = 0: ^i^^V^O — ^i^^O'^i'^iViCi^ 3; 7^ 0, y 7^ 
with the initial data T(?/, ?/, 4) = 1- At the point x = we have the jumping condition 

y, = ^3-^0 (0?^(o-, y, 6, y o, (3.23) 

where Fq{C) is the crucial gluing operator (2.5) at the defect point taking naturally the form 

Fo°(0 = ^^"\0-)M(^,a)J^(0+), where J^(0±) = exp(^^5^^^). (3.24) 

Similar to the SG case we gauge transform T T{x,y,£,) as in (3.8) and represent it as in (3.9), 
where W solves a Riccati type equation (3.11) and Z is found explicitly in terms of W as in (3.10). 
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For finding the conserved quantities C^, though we use again the same expansion at ^ — oo: 
Z{x,y,$,) = ^^^^0-3 + iJ2'?^=i ^nix,y)^~'^, the elements of the diagonal matrices -Z„(+oo,— 00) = 
Zn{+cc, 0+) + Zn{0+, —00) = Zn{+cc, 0+) + jDn + Zn{0—, — oo), have now the contribution from the 
defect point: —iDn = Z„(0+, — oo) — Z„(0— , — oo). Therefore the general form of the set of conserved 
quantities may be given by = C+ + C^ + C~, n = 1, 2, . . . , where = trace{a^Zn{+oo, — oo)), 
with 

C+ = trace(c73Z„(+oo,0+)),C~ = irace(c73Z„(0-, — oo)), C° = —itrace{uzDn). (3.25) 
Following therefore the above approach [7] with our extension, in place of (3.14) we arrive at 

Z„(+oo, -oo) = - / a2{e''''''Wn-i{x)-Wn+i{x))dx+^Dn+- / a2(e'""3W^„-i(x)-W^„+i(x))dx 
4 J-oo I 4 Jo 

(3.26) 

where Wn = —w^a^ + Wn(J- are the known solution of the Riccati equation (3.16). 

For deriving the defect contribution L>„, n = 1, 2, . . . explicitly, introduce the limiting monodromy 
matrix 

T_(x,e)= lim T{x,y,i)E_{y,^), ^heve E_{x,i) = e^-^'^^ E{x,i), E{x,i) = ^. M ei^(€-«">3a^, 

y-*-oo ^J2 \ I I J 

(3.27) 

with limy_>_(X) u{x,t) = 2Tm Using 0,(—oo) = 1 and the jumping condition (3.23) we get 

f_(0+,0 = Fo(Or_(0-,0, where r_(x,0 = n-\x)T^{x) (3.28) 

and 

= |±:^, where iJ = af ° ^"M , A = 6%, = eiiuHo)±u- m ^3 29) 
^ — la \ —A / 

Using (3.9) and (3.27) it follows from (3.28) that (1 + W{0+,^)) exp{Z{0+, -00,^)) = Fq{^){1 + 
VF(0— , ^)) exp(Z(0— , — 00, ^)) or readjusting, 

(1 + T^(0+,0)expL>(0 = FoiOil + W{0-,0) (3.30) 

where D(^) = Z(0+, — 00,^) — Z(0— , — 00,^) is responsible for generating the addition to the conserved 
quantities due to the defect. 

Note that in equation (3.30) the two unknown quantities D{^) and M^(0+, should be determined 
through two other known quantities Fq{S^) and H^(0— ,^), where i^o(0 is given explicitly as (3.29) 
and VF(0— ,^) is a solution of the known Riccati type equation for the SG model for the half-line 
X G (—00, 0). For solving -D(^) we consider expansion for large values of ^ — ^ 00: 

00 00 00 

expL>(0 = exp(5^ L>„r") = 1 + E ^^r", and 1 + W{x,0 = E ^n(x)r" (3-31) 

n=l n=l n=0 

and similarly for vanishing values of ^ — > 0: 

00 00 00 

expD(0 = exp(X^D_„n = l + E^-n^"' ^nd 1 + W(x, C) = E (3-32) 

n=l n=l n=0 
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Let us evaluate first the case with large ^ which yields from (3.30) using (3.31) the equation 



(E ^n(o+)r")(e - io){i + Yl ^ncn = + i?) e ^n(o-)r", (3.33) 

n=0 n=l n=0 

where VFo(2^) = l+ici- Gathering coefficients before different powers of in the matrix equation (3.33) 
one gets a recurrent procedure for solving the off-diagonal matrices 1¥„(0+) and diagonal matrices Dn 
from the knowledge of Wn{0—) and H: 



Wo(0+) = Wo{0-), (3.34) 

Wi{0+) - iaWo{0+) + Wo{0+)Di = Wi{0-) + iJWo(O-), (3.35) 
^2(0+) - iaWi{0+) + iWi{0+) - mWo(0+))£>i + Wo(0+)£>2 

= ^2(0-) + HWi{0-), (3.36) 



and so on. For instance, the first nontrivial result is obtained from (3.35): Di = Di = ia + iHa\ 
yielding 

= -i traceiasDi) = a{A + A'^) = 2acos W + ^ W) ^ (3 37) 

as a contribution of the defect point to the conserved quantity Cf. 

For finding next D2 use (3.36) rewriting it as ^2(0+) - iaWi{0+) + (l^i(0+) - ia{l + iai))Di + 
(1 + ia\)D2 = W2(0— ) + HWi{Q—). By taking the diagonal part of this matrix equation one gets 
D2 = HWi + iaDi, which using the relation D2 = D2 + \D\ yields 

D2 = HWi{Q-) + iaDi - ^Df, (3.38) 

where W"i(0— ) is obtained by solving the Riccati equation as wi{x) = —i{p~{x) + u~{x)). Therefore 
C2 = —i trace{azD2) is the contribution of the defect point to the conserved quantity . In this 
recurrent way we can find systematically the contribution of the defect point at x = to all higher 
conserved quantities for this integrable DSG model. Note that one can also explicitly determine from 
the above equations 

Wo(0+) = Wo(0-) = l + mi 

Wi{Q+) = Wi(O-) - aai - iaiDi + H = Wi{0-) + aiHai + H (3.39) 

etc. showing the effect of the defect on the monodromy matrix across the defect point. 

Now we switch over to the complementary case ^ — >^ and look for the conserved quantities 
„ = trace ((73 Z_„( +00, —00)) through the expansion 

Z{x,y,0 = -^^I^^3 + i E Z.n{x,y)C- (3.40) 

n=l 

We have to perform now similar expansion in the positive powers of ^ in all the above formulas 
noticing the crucial symmetry of the monodromy matrix [7] T_ ( p) = T-(x,^;u,p) , which 

is obvious from the symmetry of the SG Lax operator (2.6). It is crucial to note however that the 
unknown part of the monodromy matrix W_„(0+) across the defect point, as evident from (3.39), 
depends on both it+, u~ and obviously the above symmetry is lost. Therefore we use this symmetry 
only for u~ — ^ —u~ (with same u^) without changing the field u'^ (but with — ^ —u^, to preserve 
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the canonical structure), and expect the consistent solution of the jump condition. Therefore in place 
of (3.28) we get the condition 

r_(0+,0 =^o(e)r-(0-,O, where T_{0-,C;u- ,p-) = f_{x, ~; -u- ,p-) (3.41) 

for the known solution of the Riccati equation and 

= where = i ^ i = = exp(^(«+(0) - «-(0))), (3.42) 

Note that in (3.42) we have made the transformation ^ — |,a and u~{0) —u~{0), 

preserving p~{0) P~(0) and u^(0) u~^{0), which demands also Pq —pg for ensuring the 
corresponding quantum defect matrix Fq (5.75) to be a solution of the QYBE at the discrete level. 
This however does not affect (3.42) obtained in the continuum. 

Considering the above we obtain the corresponding matrix equations 

W.n{0+)C)h + + E ^-nD = (7 - ^) E W-n{0-)C- (3-43) 

n=0 ' n=l ^ n=0 

Arguing in a similar way we get finally the required solutions 

D-i = -ii- + Hai), (3.44) 
a 

D-2 = -HW-i{0-) - -D-i - loll, (3.45) 

where W-\ is obtained from the corresponding Riccati equation through the solution w-\ = —i(j)^{x)— 
u~{x)). Note that the contribution of the defect point to the conserved quantity C^i is 

(7°! = -itriasD^i) = --{A + A-') = -- cos (3.46) 

Qi (X Z 

while to C'L2 is — —itrace{a3D-2)- Therefore we can derive the general form for conserved 
quantities by using the simple symmetry 

C^n = (-l)"C±(p±,-^±), = (-l)-D„(-i ^z+(0),-tx-(0)), 

from those obtained in (3.26). 

Using the conserved quantities derived above we can extend now the expressions (3.21) for the 
momentum and the Hamiltonian of the SG model to include the extra contributions due to the defect 
point at X = 0: 

P(-/) = r Piu-)dx + r Piu-)dx - 2acos + 2a- cos (3.47) 

J-oo Jo 2 2 

and 

i^(-/)=r Hiu-)dx-, rg(.+)dx-(2acos"''^"^:""("^ + 2a-^ COS (3.48) 
J-oo Jo 2 2 
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where the momentum and Hamiltonian densities P{u), H{u) are given by their standard expression 
(3.21). 

To convince ourselves that (3.47, 3.48) are indeed conserved, we check it by direct calculation. 
For this we may use an identity Dt{uxUt) = Dx{^{uf + u^) + cosn), which follows easily from the 
SG equation uu — Uxx = sinu. Therefore noting that u^{x) together with their derivatives vanishes 
respectively at a; = ±oo, we get 

L>,p(<^e/) ^ 1 ((^-)2 + (^-)2) + cos u- ^{{utf - (utf) - COS u+ + (tx+ + u^p - (n+ - Ut)q) , 

\^ ^ y \x=o 



1^ 

(3.49) 



where 



. u'^ + u _i . — u 



p = asin , q = a sin . (3.50) 



Using now the Backlund gluing condition at x = 

'^t = tJ't + P + Ut=u~+p-q (3.51) 

and consequently 

(utf = {uif + + 2uip + 2utq + 2pq, 

{u+f = {u-f+p'+q^ + 2u-p-2uxq-2pq (3.52) 

we can substitute (^i^)^, {uf)"^, u+ and uf through their expressions above and apply the identity 
cosi("'" — cosn~ = —2pq to derive from (3.49) DtP^'^^f'> = 0. 

Turning now to H in (3.21) we use another identity Dt[H[u)) = Dx{uxUt), which follows again 
from the SG equation, and we show similarly that H^'^'^-^^ is also a conserved quantity. Indeed we get 

DtHi'^'^f^ = (utu- - utut + [ut + ui)p + {ut - ut)q) ^^^^ (3.53) 

where p and q are as defined in (3.50). Using again the BT (3.51) we can rewrite the first part of 
(3.53) as u^{uf — {p — q)) — ^^(^^ + (P + ?))) which clearly cancels with its second part to give zero, 
proving H^^^f^ to be a conserved quantity. 



4 Soliton solution in DSG with its possible creation &: annihilation 

We now find the relation between the scattering matrices linked to two ±-regions and the intriguing 
contribution of the defect point in creation or annihilation of the soliton. At the same time using the 
BT (2.4) unfrozen at all points as explained above we can find soliton solutions showing explicitly 
their creation, annihilation or preservation with a phase shift. 

To clarify the procedure we introduce some definitions refining that of (3.27), where wc denote 
T^^) to indicate monodromy matrix belonging to the fields u^, respectively. Remind that the fields 
have the space asymptotics: 

2'Km± for x +oo, (4.54) 
— > 27rn± for x — ^ — oo (4.55) 

which provide the following asymptotics for T^"^^ 
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t[^\x,0 ^ e^''^''^E{x,0 for x ^ -oo (4.56) 

and similarly 

^(x,0 ^ e^^3"'±E(ar,0 for x ^ +00. (4.57) 
We further relate the matrices involved using the bridging condition as 

T^^\x,y,0 = F''{x,OT^-\x,y,^)C{y,0, (4.58) 

with F'^{x, as in (2.5) and a matrix- valued function C{y,^) which does not depend on x but depends 
on y and ^. Using the chain of relations (4.54-4.57), we can relate the monodromy matrices in the it 
region as 

Tt\x,0 = F\x,OTi-\x,OF±\ where F± = diag{^ + ia±,^- ia±) (4.59) 

with 

a+ = o(-l)'"++"^-, a_ = a(-l)"++"-. (4.60) 

To get these relations one has to compare asymptotics of the functions at the infinities, choosing 
C(x,^) through F^{x,^) matrix. 

Therefore from the definition of the scattering matrix S'^^^^) = {t!^\x, ^))~^T^\x, we relate 
them as 

S^+\0 = F+S^~\OF='- (4.61) 
Now from (4.61) we get finally the relations between the scattering data 

s+ = {a+(0,?>"^(C),Ci",C2",-Cn+;7i^,72",-7X} and s" = {a-(0,i>"(0>Cr>C2"> -C_;7r,7^, -7^_} 
as 

etc., where a± as defined in (4.60) involve asymptotic im±,n±) for the fields at space-infinities and 
the defect intensity a. 

There can be three distinct possibilities [9] : 

1) a+ = (— l)™'+~'"™"a < 0, a_ = (— l)"+"'"""a > 0, when soliton number increases by 1 : A^_|- = 
N_ -|- 1 (a soliton with ^jv+ = is created by the defect). Wc have = 7" = for 
j = 1,2,...A^_, the set has an extra eigenvalue compared with S~ , and a"'"(^) = a~{$,)^^, 
6+(0 = 6-(0- 

2) a+ > 0, a_ < 0, when AT^ = iV_ - 1 and = 7^" = 7+ for j = 1, 2, ...n+, the set S" has 
an extra eigenvalue compared with S'^ and a~(^) = a'^i^^^, b'^iO = 

3) m+ -I- 171- = n+ + n_ ( mod 2), when N- = N+. The sets S'^, S~ have the same number of 
eigenvalues and ^7 = 7+ = |3^77 for i = 1, 2, ...iV+, a+(0 = a'iO, b+{0 = ?>-(e)f±|. 

In cases 1) and 2) there exists some extra defect soliton with a very special behavior. Consider 
the case 1). If this soliton moves to the right and originally is located on the left half-line x < then 
it will appear as 2-soliton after defect i.e. a soliton will be created for x > (eee Fig la). For u~ as 
1-kink solution the boundary condition (BC) gives n_ = 0,m_ = 1, while for 'u+ as 2-kink solution 
it corresponds to n+ = 0, m+ = 2, fullfilling the required condition that ri- + n+ = even, while 
m_ -|- m+ = 3, odd. In case 2) we have a similar but opposite situation. A 2-soliton moving to the 
right from x < will be converted into 1-soliton after the defect i.e. a soliton can be annihilated by 
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Figure 1: Soliton solutions (sin ' ) for DSG with a defect at x = showing a) creation, b) annihi- 
lation and c) preservation with phase shift of soliton by the defect point. 

the defect point (see Fig lb). Here for as 2-kink solution the BC can give n_ = — l,m_ = 1, while 
for u'^ as 1-kink solution it yields n+ = 0, = 1, having the required condition n_ + n-(_ = — 1 odd 
, while m_ + m+ = 2, even. 

We can derive such exact soliton solutions explicitly from the BT (2.4). For example inserting 
for u~ , 1-kink solution in the Hirota form: u~ = — 2iln j^, f± = 1 ± f, f = gl^ox+kit+if>o ^ where 
/cq = cosh0, ki = sinh0 we can extract from the BT a 2-kink solution for u'^ in the form u'^ = 
-2ilnj^, /± = 1 ± (/i + /a) + s(i(0i - 02))/i/2, fa = e'=o"'^+^^"' with scattering amplitude 
s{t) = tanh'^t, and with certain relations connecting the parameters 9, 6a and the defect parameter a. 
For A2 = — AJ = r/e*^, one gets the kink-antikink bound state (breather solution). 

In case 3) there is no creation/annihilation of soliton by the defect. In this case soliton passing 

through the defect will suffer a phase shift of cf) ^ = log ^p^, since the parameters are changed. As 

also shown in [2, 6], if we insert 1-kink u~ = — 2iln j^, /± = 1 it /i, /i = g^ox+kit+tpi g'p ^^^^ 

again have 1-kink solution for u+ = —2iln y-, /± = 1 it /a, /a = e^°^^^^^~^'^^ , with a phase shift given 
by e<^i-<^2 = _|^hd^|^ where ko = 2 cosh 0, ki = 2sinh^, a + i = 2 cosh a - ^ = 2sinhd (see 
Fig Ic). Note that the BC for the kink solutions corresponds to n_ = n_|_ = 0,m_ = m+ = 1, giving 
m_ -|- rn-4_ = n_ + n+ (mod 2), as predicted above. 

5 Classical and Quantum integrability of DSG through Yang-Baxter 
equation 

A semiclassical treatment of the DSG model through factorizable S-matrix together with some possible 
quantum features are presented in [6]. However for establishing the exact classical and quantum 
integrability, it is necessary to show the validity of the Yang-Baxter equation for this model both at 
the classical and the quantum level. Our aim is to carry out this program by finding the associated 
quantum and classical i?-matrix and the lattice regularized Lax operators for this system including 
the defect point, as exact solutions of the YBE. Subsequently we formulate the algebraic Bethe ansatz 
for the quantum DSG model. Our strategy is to follow closely the approach of the standard quantum 
SG model [11] in combination with the ancestor model scheme of [10]. 

5.1 Exact quantum integrability of lattice DSG model 

We try to construct first an exact lattice regularized version of the quantum DSG model through a 
discrete monodromy matrix 

m = T^+(Oi^o )T^-(0 (5.63) 
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where 

T^+(0 = U+i^,u+) ■ ■ ■ U+{^,ut), T^-(0 = UZ,{^,uZ^) ■ ■ ■ U:^{^,uZr,) (5.64) 

with Up{^,u^), j = ±1, ... lb TV being the discrete quantum Lax operator of the lattice SG model 
defined along both sides of the defect, while Fq{^, Uq,Uq) is the quantum Lax operator at the defect 
point j = 0. Recall that [11] for quantum integrability the monodromy matrix of the system (5.63) 
must satisfy the global version of the quantum YBE (QYBE) 

R{^,v)T{0»T{n) = {T{r,)^I){TiO^I)R{^,v), (5-65) 

which taking trace from both the sides yields evidently the relation [T(^),r(?7)] = 0, where r(^) = 
traceT{^) = J2n^n^"^, giving finally the quantum integrability condition through the commuting set 
of conserved quantum operators as : [Cn,Cm] = 0. If we ensure ultralocality condition, i.e. all 
constituent Lax operators Up, j G [ibl,ibiV] and Fq, in (5.63) situated at different lattice sites 
mutually commute, then it follows from (5.65), that each of these local Lax operators Up, Fq must 
also satisfy exactly the local QYBE: 

R{^,ri)Lj{0 ® Ljiv) = iHv) ^ I)iHO ® - V) (5.66) 

with the same quantum i?-matrix, where Lj = Up at j € [±1, ibA^] and Lq = Fq at the defect 
point. The transition from local QYBE (5.66) to the global one (5.65) is a reflection of the coproduct 
structure of the Hopf algebra property of the underlying quantum algebra [12]. 

It is known that the standard SG model satisfies (5.66) with trigonometric i?*^*^-matrix [11] 

R\\ = Rll = a(A, a), R\\ = RH = 6(A), R^ = R^ = c{a) 
where a(A, a) = sin(A + a), 6(A) = sin A, c{a) = sin a. (5.67) 

Our task therefore is to find the discrete Lax operators Up and Fq, which in one hand would satisfy 
the QYBE (5.66) with i?*''*^-matrix (5.67) and on the other hand would recover the DSG field model, 
we have started with. That is the construction of the exact Lax operator solutions of the QYBE: Up 
and Fq would be such that they would reduce in the continuum limit (lattice constant A ^ 0) to 
the field Lax operator of the SG model U{x,$,) (2.6) and the BT matrix Fq (2.5), respectively. For 
such a construction we turn to the general scheme of [10] , where it was shown that any Lax operator 
satisfying the QYBE (5.66) with i?*'''^-matrix (5.67) can be obtained as a particular realization of the 
ancestor Lax operator 

Ttri9(c^_( e4^^e-^'+r'4"^e--' 2sina4-) \ 

^'^'^^^^^ - Iv 2sina4+) CS^^e-^-^ + ) ' ^ ' 

with the quantum spin operators generating a generalized quantum algebra 

[sW, 4-)] = (mW sm{2as^) - iM^-^ cos{2as^)) [s^ s^] = ±s^p\ [M^, •] = 0. (5.69) 

Here the deforming operators M^^) = {c^^C2 ^ ±c^ ^^2^^) are expressed through a = 1,2, which 
are mutually commuting and central (superscripts (±) here are obviously different from the field labels 
zb used in the DSG). 
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We intend to construct the quantum Lax operators Up and Fq as a solution of QYBE in a unified 

way from the same (5.68). Note that a reduction as = ^iA,a = 1,2, takes (5.69) to suq{2) 
algebra 

[4+),4-)] = 2A^^^^^, [^4±)]=±4±), , = e- (5.70) 
siri cx 

which can be realized in canonical variables [u^,p^] = iSjk as 

s' = ^, s(q'^Hu^,p^) = e-'^^^g{u^,A), 4") = (4+))t, (5.71) 

where c/(n. A) = (1 + cos Q;(n + 1)) ^ 

Inserting the generators of suq{2) realized as (5.71) in (5.68) we recover the quantum Lax operators 

jj± ^ I sina(A + ^) smas^f\uf ,pp \ ^^^iaX ^572) 

\ sin as'^\uf ,pf) sina(A-^) / 

for the fields on the lattice in conformity with the SG model [13]. 

For constructing the discrete BT operator Fq{^, u^,p^) again from (5.68), we choose the reduction 
ci^^ = 1, ci ^ = 0, a = 1,2, giving M^^^ = and reducing (5.69) to a simpler algebra 

[4+), 4-)] = 0, e^"^%(±) = e±^"4^^e'"''- (^■'''3) 

Fortunately, we can find a consistent realization of the algebraic relations involving both the canonical 
fields {w^ ,p^) in the form 



e^"^' = e_, 4+) 



(4"^)^ = ae+ Pr\ where e± = e^f ("0 ±«o ), p_ = e^^CfJ-Po ). (5.74) 



With commutation relations [e+,P_] = 0, e_P_ = e~'"P_e_, it is instructive to check that the 
generators (5.74) satisfy the algebra (5.73). Therefore the Lax operator (5.68) reduces finally to the 
explicit form 

i^o(^,4,Po ) = Pt'F^it4)P:^'"\ F^i^,u+,u-) = e-x-3«-M(e,a)eX-3«+, (5.75) 

connecting remarkably to the BT operator Fq{^,u^) (2.5) for the DSG model. Note that both the 
above discrete Lax operators obtained as realizations of the quantum intcgrable L-opcrator (5.68), by 
construction must satisfy the QYBE exactly with the i?*'"*^-matrix. Consequently, (5.63) with (5.64) 
represent a quantum integrable discrete DSG model. 



5.2 Exact classical integrability 

As the quantum Lax operator should satisfy the QYBE for its quantum integrability, the corresponding 
classical Lax operator consequently should satisfy its classical analog the classical YBE (CYBE) [14] 

{Lj{0,^Lj{n)}pB = 6jk[r{^,r]),Lj{0 ® L^iv)]. (5.76) 

Note that at the classical limit h ^ 0, the quantum commutator should reduce to the Poisson bracket 
(PB) and the operator elements of the quantum Lax operator would become just functions, with the 
form of the Lax operators (5.72, 5.75) remaining the same. Only the normal ordering needed in the 
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quantum case should be ignored now. However since the i?-matrix provides the structure constant for 
the commutation relations, for transition to the classical limit we have to scale the parameter a in the 
i?-matrix (5.67) as a ha. This therefore defines the classical limit as a ^ in all the elements of 
the quantum i?-matrix reducing it to the classical r-matix: = I + ar{X) + o{a), r(A) = 

^j^daR\a=o as 

rll = ri = ao(A), r^^ = rjj = 0, r^l = rfl = co(A) 

where ao(A) = cot (A), co(A) = (5-77) 

sm A 

while the QYBE (5.66) reduces to the CYBE (5.76). The underlying algebraic relations are turned 
into PB relations at this classical limit, which is achieved by taking a ^ in (5.69), only in the terms 
that come from the i2-matrix. As a result the classical limit of (5.69) is given by 

{s'^q^\4~^}pB = (m^"^^ sin(2as=^) - iM^'^ cos(2as^)) {e*"^', s^Ipb = ±^e^"^%W, (5.78) 

Therefore similar to the quantum case, any realization of the classical limit of the L-operator (5.68), 
with its elements satisfying the PB relations (5.78) must be a solution of the CYBE, by construction. 
Consequently, since (5.72) and Fq ( 5.75), as shown above, are indeed realizations of (5.68), 
their classical limit must satisfy the CYBE (5.76) exactly, proving the classical integrability of the 
DSG model. Alternatively one can check this statement through direct verification. For example at 
M(+) = 2A2,M(-) = 0, i.e. at the sug{2) limit the PB relations (5.78) reduce to 

H-'\s(-^}ps = 2A^'^^^, {e-\s^^^}pB = ±^e-^%W, (5.79) 
^ ^ sin a ^ ^ 

which is evidently satisfied by the realization of the generators in the canonical fields {u^,Pk}pB = ^jk 
as ^ 

s^ = —, 4±)(u±,p±) = -^e^^'P^ (1 + A\os au)^ (5.80) 

2 ^ sin a 

using the definition {f,g}pB = ~ i^lf- * similar way we can check the validity of the required 

PB relations for the elements of (5.75), obtained from (5.78) at M(=^) = 0: 

{4+\ 4"^}PB = 0, {e*"^', 4±)}PB = ±ie-^%(±). (5.81) 

Using the classical analog of (5.74) expressed in the same form without normal ordering: e^"'^ = 
e_, Sq^^ = ae+ -P?^, and the PB relations like {e+,P_}pB = 0, {e_,PT^} = ^ie^,P^^ we verify 
the relations (5.81), which guarantees that Fq satisfies the CYBE (5.76). The solution of the local 
CYBE (5.76) leads also to that for the global CYBE for T(0 = Uj Lj{0- 

{T{0,mv)}pB = Sjk[ri^,v),TiO^T{n)], (5.82) 

from where in exact analogy with the SG model [7] one can extract the action-angle variable corre- 
sponding to the continuum as well as the discrete spectrum, proving the complete classical integrability 
of the DSG model. This also reveals an intriguing fact that the local differences between the SG and 
the DSG models seem to become irrelevant at the global level (5.82). 

It is remarkable that at the continuum limit A ^ 0, when the discrete Lax operator Fq Fq 
recovers the BT matrix (2.5), the dependence of the canonical momentum P_ drops out completely. 
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creating a paradoxical situation that would result a trivial PB for the F° and naturally not satisfying 
the CYBE (5.76). Therefore for showing the validity of the CYBE for the field model one has to be 
careful and should take first the PB in the corresponding exact discretized model and then perform 
the continuum limit. Similar is true also for the Lax operators U^. 

Our another important result is a discretized BT relation that we can derive connecting the discrete 
sine-Gordon Lax operators (5-72) in the form of a discrete gauge transformation 

C/+(e,u+) = F/+i(e,u±+i,p±,i)C/r(e,uT)(F/(e,ix±,p±))-S (5.83) 

through the gauge matrix Fj'{^, u^,pf) (5.75). Note that we assume this discretized BT to be valid for 
arbitrary j. At the continuum limit the discrete gauge transformation (5.83) should recover the known 
gauge relation (2.4) connecting the field Lax operators U^{u^). Similarly comparing the elements of 
the matrix relation (5.83) we can obtain the bridging relation between the discrete variables and 
uj of the discretized DSG model at any site j, involving also the variables 77,^^,p^,p^_^^. This discrete 
BT relation would yield at the continuum limit A — >■ 0, the standard BT relation between the fields 
u^{x) (2.2). 

5.3 Algebraic Bethe ansatz 

Following the formulation of quantum SG model [11] we can apply the Algebraic Bethe ansatz 
method to the lattice regularized quantum DSG constructed above and solve in principle its eigen- 
value problem exactly. Recall that the aim of the algebraic Bethe ansatz is to solve exactly the 
eigenvalue problem of r(^) = traceT{S,), r(^) = YljLj{S^), generating all conserved operators in- 
cluding the Hamiltonian, in the form: r(,^)|n >= A(^)|?7, >, with the eigenstates \n > defined as 
\n >= 1^1, . . . >= n" ^(6)|0 >• TuiO = -8(0 acts as creation operator, while T2i{(,) = C{^) as 
destruction operator annihilating the pseudovacuum: C{$,)\0 >= 0. A crucial step in this formalism is 
to construct the pseudovacuum state |0 >, which is achieved for the SG model by combining the ac- 
tions of consecutive pair of Lax operators: UjUj+i\0 > [11]. Repeating this procedure we can construct 

±— (2) 

the pseudovacuum |0 nj=±i \^j > foi" the quantum DSG model, yielding C=''(A)|0 0, for 
all sites except the defect point at j = 0. 

However the single defect point would play a nontrivial role, since after crossing this point, say 
from the left the pseudovacuum property: 

.-(0,0 >-^(-<«l° 

would be lost due to nontriangular matrix form of FqIQq >. Instead of annihilating the local vacuum, 
as needed by the lower left corner operator element of the matrix Fq\Qq >, the defect at site j = 
would turn it to a state 0\i^o >, O = — ae**^^*-^" ^^i' +"0 ))^ creating at the same time its conjugate 
state -O^no > at the upper right corner. This is expected to lead to the creation/annihilation of 
quantum states by the defect point similar to that with classical solitons as we have observed exploiting 
the BT in sect. 4. Perhaps one should explore the possibility of using a quantum extension of the 
BT [15] to mimic the successful classical approach, using the relation like (5.83) at the quantum level. 
These tricky points however need careful and separate analysis and should be dealt with elsewhere. 
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5.4 Continuum limit 



It is crucial to check that the classical and quantum integrable discrete DSG model we constructed 
and solved above would yield the same DSG field model we have started with, at the continuum limit: 
A — 0. Note that at this limit the canonical variables go to canonical fields: u'^{x),p^ 
Ap^(x), with [u^ (x) , (y)] = id{x — y). Therefore for extracting the limit we have to scale giving 

e"^^ ~ 1 + iAp^{x), in both the discrete Lax operators and Fq. 

Observing further that rotated operator (JiUj'''^{X) is also a solution of QYBE due to the symmetry 
of the i?-matrix: [i?*"^, (Jq (g) cja] = 0, a = 1, 2, 3 , we expand the Lax operator (5.72) in powers of A 
at all sites j (except at the defect point) as 

a^Uf = (1 + / ' U^{x')dx') + 0(A2) = (1 + AU^{x)) + 0(A2), (5.85) 

with the field operator U^{x) recovering exactly the Lax operator (2.6) of the SG field model. Note 

that we have put the deforming parameter a = 1 in all expressions related to the continuum model, 
for simplicity. Thus we recover at the continuum limit, the DSG field model at all points except at 
a; = 0. 

Performing the same continuum limit at the defect point j = , we get on the other hand from 
(5.75) : 

Fi{i,u+,p+-u^,p^)^F^{i,u+{Q),u-{Q)) + iAF', f' = {p+ -p^)[as,F^{^,u+{0),u-{0))], (5.86) 

clearly giving F° at A — >^ 0, i.e. it recovers the same BT matrix (2.5) at the defect point, meeting the 
essential requirement. Therefore collecting all nontrivial terms we get finally the continuum limit of 
the lattice regularized model (5.63) as 

T(A)= (e/o^°°^^(^'"')''"') FO(A,u+(0),n-(0)) (^e-/'-°o (5.87) 



yielding the original DSG field model. 

Finally in the continuum limit using (5.85) and the expansion like F^_^i — > F°+A(F^+iF'), u^^-^ 
u^ + Au^, pf^i A(p^ + Ap^) etc. we can show directly that the discrete BT relation (5.83) goes to 
the field BT as gauge-transformation between U^{u'^), while the relation between its matrix elements 
connecting recovers the bridging relation (2.2). 

We stress again that for proving the classical and quantum integrability of the field model (5.87) 
one has to go to its proper lattice regularized version (5.63) and check the validity of the GYBE and 
QYBE before taking back the continuum limit. 



6 Concluding remarks 

We have proved here the classical and quantum integrability of the sine-Gordon model with a defect 
by finding the exact solution of the related quantum as well as classical Yang-Baxter equation through 
integrable discretization of the model. This enables one to solve the quantum eigenvalue problem by 
algebraic Bethe ansatz and the classical model for the action-angle variables. 

Combining Matrix Backlund transformation with the matrix Riccati equation approach we have 
extended the existing formalism and found all higher conserved quantities for the defect sine-Gordon 



16 



model in a systematic way, with explicit forms of the defect contribution. In the simplest case this gives 
the momentum and Hamiltonian of the DSG model found earlier. Deviating from the earlier studies we 
have used the unfrozen BT, which can produce intriguing effect like creation or annihilation of solitons 
by the defect point, apart from the preservation of soliton with phase shift as predicted earlier. We have 
also found the constraints showing exactly how the creation/annihilation or preservation of solitons 
depend on the boundary conditions of the field, in the framework of the defect sine-Gordon model. 
This result should be of crucial importance in the possible experimental detection of such unusual 
events. A pertinent question arises here regarding the obvious violation of topological charge in this 
SO model with a defect due to possible nonconservation of soliton number. It should be noted however 
that the topological charge arises in the SO model as a degree of mapping from S^, while a 

defect in the coordinate-axis or a discontinuity point (like a puncture in the sphere) as occurs in the 
DSG, can not be mapped into a smooth sphere or S^, violating thus the concept of the topological 
charge itself. Therefore in the DSG the solitons seem to be no longer topological and hence their 
number may change. The formation of nontopological semi-fiuxons [4] therefore may also be possible 
to explain in the framework of the DSG model [16]. 

After completion of our work, very recently an important paper has appeared in the arXiv [17], 
where systematic studies were made using the Lax pair formalism of integrable systems in the line 
of the present investigation, and consequently infinite set of conserved quantities were obtained for 
a whole class of defect models e.g. SG, NLS, KdV, mKdV, Liouville, DNLS etc. , belonging to the 
AKNS and the KN spectral problems and having a defect at a single point, which in principle could be 
extended to multiple defect [17]. Similar to our approach the conserved quantities are found through 
Riccati equation, though in contrast to the Matrix Riccati equation used here this equation is a scalar 
one. In the same work the importance of establishing the complete integrability of the defect models, 
at the classical and the quantum level, through the Yang-Baxter equation as well as the necessity of 
discretization of the model for achieving this goal are emphasized. Interestingly both these agenda in 
the wish-list were addressed and solved rigorously in the present paper. 

Extension of the exact result of the quantum and classical DSG model through the Yang-Baxter 
equation presented here, to other models like NLS, DNLS etc. treated in [17] would be an interesting 
problem. 
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